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Motivated by recent experiments on toroidal Bose-Einstein condensates in all-optical traps with 
tunable weak links, we study the one-dimensional Bose-Hubbard model on a ring-shaped lattice 
with a small region of weak hopping integrals using quantum Monte Carlo simulations. Besides the 
usual Mott insulating and superfluid phases, we find a novel phase which is compressible but non 
superfluid with a local Mott region. This 'local Mott' phase extends in a large region of the phase 
diagram. These results suggest that the insulating and conducting phases can be tuned by a local 
parameter which may provide a new insight to the design of atomtronic devices. 
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Introduction. Cold atom experiments utilizing an op- 
tical lattice provide an excellent testbed for quantum 
many body problems which were previously inaccessi- 
ble in conventional materials. A remarkable achievement 
is the realization of the Bose-Hubbard (BH) model us- 
ing ultracold atoms on optical lattices^ with the addi- 
tion of a confining potential that results in the "wedding 
cake" structure^. Over the last two decades, a consider- 
able amount of work has been devoted to understand the 
ground state phase diagram of the BH model and its vari- 
ants— ~— . In general, the model contains a superfluid (SF) 
phase at incommensurate fillings and a Mott insulating 
(MI) phase at commensurate fillings and strong coupling. 
The SF phase is gapless, whereas the MI phase is char- 
acterized by the existence of an energy gap for creating 
a particle-hole pair. As the density is changed or the in- 
teraction strength is varied, the BH model can be tuned 
from the MI to the SF. Tuning between insulating and 
conducting phases by controlling the external parameters 
provides a tantalizing opportunity of creating analogs to 
electronic devices and circuits by using ultra cold atoms 
in optical lattices, which have been recently defined as 
'atomtronics— The conventional electronic system is 
based on the electron charge, whereas the atomtronic 
system can use neutral atoms which are either bosons 
or fermions, moreover the optical lattice is better con- 
trolled. Based on this unique property, it has been sug- 
gested that these atomtronic systems may be useful in 
quantum computing™.. Some theoretical models have al- 
ready been proposed for atomtronic devices such as bat- 
teries, wires, diodes, and transistors^— . 

A recent advance on optical lattices is the realization 
of confining potentials with toroidal shapes by using the 
intersection of two different red-detuned laser beams^^. 
The versatility of this technique allows the creation of 
ring-shaped lattices by superimposing an optical lattice 
on a toroidal confining potential, which is a realization of 
a quasi one-dimensional lattice with periodic boundary 
conditions. Remarkably, it is possible to control the local 
hopping parameter in a region of the ring by applying 



a magnetic field and an additional laser bearn^. This 
opens up the new possibility that the different phases 
in a boson system not only can be tuned by a global pa- 
rameter, such as the coupling strength or chemical poten- 
tial, but also by a local parameter, such as the tunneling 
strength of a small region of the entire lattice. It has 
been suggested that this property can be utilized as an 
alternative realization of atomtronics«. 

In this letter, by using quantum Monte Carlo (QMC), 
we show that introducing weak links in a ring lattice 
can produce a local Mott (LM) phase in addition to the 
usual MI and SF phases present in the homogeneous BH 
model. Zero temperature local incompressible MI behav- 
ior was shown in a one dimensional system of interact- 
ing bosons in a confining potential^. Our non-confined 
model exhibits a LM phase which is gapless and non-SF, 
and a region of LM insulator which exhibits incompress- 
ible MI behavior. This is an important result which sug- 
gested that by controlling the local tunneling strength 
the system can be tuned between a SF phase and a MI 
phase thorough a non-SF LM phase. This provides the- 
oretical support that atomtronic switches can be imple- 
mented by tuning certain local parameters in a quasi one- 
dimensional system. 

Model and method. We consider a bosonic system on 
a torus-shaped lattice, where the section of the torus is 
sufficiently small compared to the primary radius so the 
physics can be reduced to a one-dimensional lattice with 
periodic boundary conditions. The Hamiltonian takes 
the form 
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where L is the number of lattice sites. The creation and 
annihilation operators a\ and a i satisfy bosonic commu- 
tation rules, [a^aj = [af,at] = 0, [a^aj] = 6^, and 

hi = a\a i is the operator that measures the number of 
bosons on site i. The parameter t is the global mag- 
nitude of the hopping integral. In this paper, we use 
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t = 1 to set the energy scale. The sum runs 0VG 

all distinct pairs of first neighboring sites i and j, anc 
Wij £ [0; 1] determines the weakness of the hopping inte 
gral between i and j. In the following we consider a sys 
tern with M consecutive weak links for which Wij = J j t 
where J S [0; t] is a control parameter, and L — M stronj 
links with m,j = 1. We restrict our study to the casi 
with 10% of weak links (M = L/10). The parameter I 
determines the strength of the on-site interaction. 

In order to solve this model, we perform exact QMC ii 
both canonical and grand-canonical ensembles by usin; 
the Stochastic Green Function algorith m 15 ! 16 with globa 
space-time updates^. In the canonical ensemble, th< 
number of particles N is a parameter and remains con 
stant during the simulation. The chemical potential p is 
measured at zero temperature by the finite energy differ- 
ence p(N) = E(N + 1) — E(N). In the grand-canonical 
ensemble, the number of particles is given by the quan- 
tum average of the operator N = J^. hi, and is controlled 
by adding to the Hamiltonian (TTJ) the term —pN where 
p is a control parameter. We use an inverse temperature 
/3 = L/t in order to capture the ground-state properties. 

Superfluid density and compressibility. For the uni- 
form system, J = t, only two phases are present: Mott 
insulator (MI) and superfluid (SF). The MI phase oc- 
curs at commensurate fillings and large onsite repulsion 
U , and is characterized by a vanishing compressibility, 
dp 



dp 



where p = N/L. The SF phase is detected 



by measuring the superfluid density, p 8 , given by the r 
sponse of the system to a phase twist of the wave functic 
at the boundaries. In our QMC simulations, it is conv 
nient to relate this superfluid density to the fluctuatioi 
of the winding number, W, via Pollock and Ceperley 
(W 2 )L 



formula—, p. 
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We have checked analytical 



and with exact diagonalization that the above formu 
remains valid for the non-uniform system where J < t. 

In the following, we show that there exists a ran; 
of parameters for the non-uniform system for which \ 
observe a vanishing superfluid density and a finite cor 
pressibility at incommensurate fillings. Fig. [1] shows tl 
superfluid density p s as a function of the chemical pote: 
tial p for L = 50 and {7 = 8. Here we use grand-canonic 
simulations for different weak link hopping J. We ce 
clearly sec that the region with vanishing superfluid den- 
sity expands over a large range of chemical potentials p 
when the strength of the weak links is lowered. 

Fig.[2]shows the density p and the superfluid density p s 
as functions of the chemical potential p, for both homo- 
geneous ( J = t) and inhomogeneous ( J = 0.2t) systems, 
for L = 50 and U = 20. We can see that a Mott plateau 
at p = 1 exists until p = 16.1 with a vanishing super- 
fluid density p s and compressibility k, for both systems. 
For p > 16.1, the density p starts to increase and the 
compressibility k is finite. As it is well known, the super- 
fluid density p s of the homogeneous system is non-zero 
as soon as the density is no longer an integer. However, 




FIG. 1: (Color online) The superfluid density p 3 as a function 
of the chemical potential p for L = 50 and U — 8 in the 
ground state. The figure shows results for different values 
of the weak hopping integrals, J = 0.2t (circles), J = OAt 
(stars), J = 0.6t (triangles) and J = 0.8t (diamond). 



for the inhomogeneous system, the superfluid density re- 
mains zero until p re 19.1. Thus there exists a finite 
range of values for the chemical potential for which the 
superfluid density is vanishing but the compressibility is 
finite. Therefore, as the chemical potential is increased, 
the inhomogeneous system undergoes a phase transition 
from a MI phase to a new phase, then to a SF phase. 
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FIG. 2: (Color online) The density p and the superfluid 
density p s as functions of the chemical potential p for the 
homogeneous system (J = t) and an inhomogeneous system 
(J = 0.2t), for L = 50 and U = 20. 

Properties of the phases. We investigate the interme- 
diate phase first by analyzing the local density of the 
lattice. The local density in the homogeneous model is 
uniform, whether the system is in the MI or SF phase. 
For the inhomogeneous model, we have phases with a 
non-uniform local density, as shown in Fig. [3^,. We in- 
sert 10% of the weak links in the middle of the lattice. In 
the MI region, the local density n% throughout the entire 
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lattice is uniform and sticks to integer values (rij = 1 for 
the first Mott lobe, n i = 2 for the second one, etc). 

When additional particles or holes are added to the 
lattice the weak link region keeps its integer density (see 
Fig. [3k). Outside the weak link region, the local density 
shows an oscillatory behavior. These two observations in- 
dicate that the additional particles do not affect the MI 
character of the weak link region until the number of ad- 
ditional particles or holes is beyond a critical density. For 
a one-dimensional system, the superfluid density or the 
winding is zero when part of the system is locally Mott. 
As a result we identify this locally integer-density region 
as a local Mott (LM) phase. The weak link provides a 
fixed boundary condition for the density profile, the ad- 
ditional particles or holes accumulate outside. Then, the 
region with J = t can be effectively described by the 
hard-core boson model with L — L wea k number of sites. 
For a one dimensional system, the hard-core boson can 
be written in terms of spinless fermions using the Jordan- 
Wigner transformation,— and the oscillation of the local 
density can then be explained by Friedel oscillations^ 
where rii ~ cos(kFXi), where kp is the Fermi wavevector 
given by the particle density. This explanation is corrob- 
orated by the numerical data which show that the cycle of 
the oscillation of the local density is approximately given 
by l/|n-1.0| for \x = 16.6 and l/|n-2.0| for fj, = 22 (see 
Fig. Hi). 

When adding more particles beyond the critical den- 
sity, the local density at the weak link shifts away from 
integer values. This suggests that the LM insulating re- 
gion is destroyed. Thus it opens the path for the flow, 
and we find that the superfluid density becomes finite 
when this happens. 

We study the dynamics of the model by evaluating 
the low energy excitation spectrum. Using the Feyn- 
man single-mode approximation, the low energy excita- 
tion spectrum Q(k) can be written as^i 



fi(fc) 



E k 
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where, 



E k = — (cos/c - 1)(* |X]KWi + 4fi a i)l*o), (3) 



l^o) is the ground state, and S(k) is the static structure 
factor. 

Fig. [3Jd displays the low energy excitation spectrum 
throughout the reciprocal lattice space. In the MI region 
it shows a gap near zero wave vector, whereas it has a 
linear dependence for the SF phase. The linear behavior 
is expected in the SF region due to the gapless Gold- 
stone mode. In the LM region the low energy spectrum 
shows a parabolic behavior, as expected for disordered 
free-particles. Since the LM does not follow a linear be- 
havior near k — 0, no signal of super- flow exists in the 
LM region. 
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FIG. 3: (Color online) The local density (top panel) and the 
excitation spectrum (bottom panel) for L = 50 and U = 20, 
in the ground state. Top panel: The local density m as a 
function of the site index i for different values of the chemical 
potential /i. Bottom panel: The low energy excitation spec- 
trum Q(k) in the three regions of the phase diagram: SF, LM 
and MI. 



Ground state phase diagram. The MI phase is charac- 
terized by an integer local density and the existence of a 
finite gap for single particle excitations. At zero temper- 
ature, the gap can be easily obtained in the canonical en- 
semble. We define the gap as A = fi + — where ^_ and 
/i+ are the minimum and maximum values of the chemi- 
cal potential for which the MI phase exists. By definition, 
M+ = E(N+1)-E(N) and n_ = E(N)-E(N-1), where 
N is the number of particles in the MI phase. The func- 
tions n~(t, J,U) and fi + (t, J 7 U) determine the bound- 
aries between the MI and LM regions. Since the total 
density remains unchanged for fj, G [/z_;/z+], the com- 
pressibility k is vanishing in the MI region. 

We determine the phase boundary between the LM 
and SF regions by using the grand canonical simulations 
and scanning over the chemical potential, as in Fig. [2] 
The critical value /i c of the chemical potential where the 
superfluid density becomes non-zero depends on the size 
of the system, L, and converges to a finite value in the 
thermodynamic limit. As the size increases, the curve 
displaying the superfluid density becomes sharper and 
sharper. Since we work with a fixed large size, L = 50, 
we define /x c by the value of the chemical potential that 
corresponds to the maximum slope for the superfluid den- 
sity curve. The curve /z c (t, J, U) determines the bound- 
ary between the LM and SF phases. 
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FIG. 4: (Color online) The ground state phase diagram ol 
the inhomogeneous system (J = 0.2t) in the (u/U, t/U) plane. 
The lines with solid circles show the first and the second Mott 
lobes, and the lines with open circles show the boundaries 
between the LM and the SF regions. 



In our simulations the density varies continuously as a 
function of the chemical potential. This suggests that the 
transitions from MI to LM and from LM to SF are contin- 
uous, as it is the case for the homogeneous model22r— . 
We show in Fig. [4] the ground state phase diagram for 
J = 0.2f in the (n/U,t/U) plane. The Mott lobes that 
are present in the homogeneous model are weakly de- 
formed by the presence of the LM phase. The phase 
boundaries near the tip of the Mott lobes are difficult to 
estimate due to the very small LM region. 

We now investigate the variation of the phase bound- 
ary between the LM and SF regions as a function of weak 
link hopping J, Fig. [5] For a fixed value of the interaction 
U, the phase boundary lifts up linearly when decreasing 
the hopping J in the weak link reducing the size of the 
SF region in the phase diagram. In the limit J = the 
curve extrapolates to /i/U = 1, and the SF region com- 
pletely disappears since the system is no longer periodic. 

Conclusion. In this study, we propose that a super- 
conducting ring with weak links might display a new 
phase which is gapless, compressible, and non-supcrfluid, 
with local Mott insulating behavior. This phase does not 
exist in the homogeneous Bosc-Hubbard model. In the 



experiment by Raman et al., a toroidal condensate is 
created with a smooth trapping potential^. If the ex- 
periment can be repeated by superimposing a lattice on 
top of the toroidal potential, our model could be directly 
studied experimentally. Our results have direct impli- 
cations for atomtronic device s 11 ' 13 . For example, if the 
chemical potential /i, the weak link hopping J, and the 
interaction U/t are tuned so that only the link is a Mott 
insulator, then a gate above the link can be used to switch 
the conductivity of the link on and off. The non-linearity 
of the switching can be tuned by adjusting the link width 
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FIG. 5: (Color online) The critical value of the chemical 
potential /j, c between the LM and SF regions, as a function of 
the weak link hoping J, for U — 8 and U = 20. The variation 
with J is quasi-linear. 



and hoping J/t. Complex circuits with highly non- linear 
behavior may be constructed by a series of such switches. 
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